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1. Introduction 


The objective of our Investigations Is to analyse the magnetic dis- 
turbance fields caused by global external current systems with particular 
enq;)has1s on Improving the understanding of the physical processes which control 
high latitude current systems. These objectives and approach have been outlined 
in detail In Triennial Progress Report #1 (December 1, 1980). 

Work towards these objectives has progressed satisfactorily, as 
outlined below. 

2 . Techniques 

2.1 Data Display 

Following processing of the data as described In Report #1, the 
data Is routinely plotted In the Universal Time (UT) format (Figure 1) as well 
as In a polar plot format (Figure 2). These plots have now been standardized 
to consist of the data In the K'D'U' coordinate system (see below). The polar 
plot of the H'D' components also has a panel showing the U' component. Both 
standard plots are labelled with EKp, Kp, Dst and the X, Y and Z components of 
the Interplanetary magnetic field (IMF) corresponding to the time of the highest 
latitude of that transpolar segment of the orbit. 

2.2 Coordinate System 

The H'O'U' coordinate system has been adopted as the standard for 
representing the Magsat residual magnetic field vectors (I.e., the measured field 
- 6SFC MG5T 06/80 model field). This system Is defined as follows: 

Let E = unit vector In the direction of the eccentric dipole axis derived from 
the value of the coefficients g^, gj and h{ in the MOST 06/^ model. 

Let 0‘ = unit vector tangent to the MGST 06/80 model field at the point of obser- 
vation. 0* is antiparallel to the field In the southern hemisphere, 
parallel to the field in the northern hemisphere. 

Let Z‘ = unit vector from the point of observation to the eccentric (I.e., offset) 
dipole origin derived from the coefficients g2, g^ and hi In the MGST 
06/80 model. 

Then H' « U' x (Z' x E); note that H' lies In the meridional plane of the eccen- 
tric dipole and is perpendicular to the model field vector. H' always 
has a northward component, but Is not usually horizontal. 

A A A A 

0' * U' X H'; D' is always eastward ^directed but is not necessarily horizontal, 
because of the declination of U*. 

It is felt that a field-aligned system Is more natural than a verti- 
cally aligned system for spacecraft data b<icause of the sensitivity of high 
altitude magnetometers to field-aligned current (FAC). 
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2.3 Data Selection 

It has become necessary to select subsets of the Magsat data set for 
particular studies. For exan^1e» to date» only hourly averages of the frtMn 
the ISEE-3 satellite have been available. It has been found In preliminary 
correlative studies of the polar cap magnetic field and the IfV By that most 
significant results are obtained during periods of high I^^ stability. The 
“stability" Is represented by a five-digit nuiimr derived from parameters from 
the ISEE-3 data tape, as follows: 

First (highest order) digit = 1 1f the sign of Bx, By and B 2 of the IMF (In GSE 
coordinates) during the IJT hour of 1nter4st Is the same as during 
the previous hour. 

» 0 otherwise. 

Second digit = IP » [lOf] 

where [] = "largest Inte^r in" 

f = fraction of milliseconds in the UT hour for which IMF 
data exists 


and IP is set to 0 if IP < 0 
IP is set to 9 if IP > 9. 

Remaining digits = [lOd]: where d = {(<B,->)2/<Bj2>}* 

for digits 3, 4, 5 and i = x, y, z (in GSE coordinates): 
where <> indicates hourly averages. 


)lishments and Work in Progress 


3.1 Data Management 


a) A data file containing the following parameters has been generated 


i) IMF Bx» By, B, (GSM coordinates); hourly averages from ISEE-3 only 
and IMF stability 

ii) Ost, and daily average Dst. The values are those provided on 
Sugiura's tape, augmented by NASA-Goddard values when the Sugiura 
data are not available 

iv) 

v) 

vi) 3- hour Kn 

vii) daily EKp 

b) A computer code (TPOLAR) has been developed to generate, from the orbital 
elements, the time, latitude and MT of the extremum latitude of each 
transpolar se^nt of orbit. The results of TPOLAR are stored on disk. 
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c) A computer code (SEARCH) hes been developed to Identify passes (trans- 
polar orbit segments) for which any of the Indices described In a) fall 
Into any desired range. It Is possible for example to list all passes 
In the southern hemisphere with By < Of 6^ > 0 and f10.7 > 200. 

d) A counter code Is being developed to determine equatorward and poleward 
boundaries of the regions of transverse perturbation vectors which are 
thought to characterize the auroral oval. The code examines the direc- 
tion and magnitude of the hi^ latitude perturbation vectors. The 
criteria for boundary specification are being refined at this time. 

3.2 Modelling 

The existing magn f field modelling code, using a refinement of 
Kisabeth's grid cell system (see Report #1), is being used to examine theoretical 
models of auroral current systems and to predict Magsat observations. The com- 
puter code in use at NRCC permits specification of the current strength with a 
resolution of 0.5° in latitude and 4° in longitude, with a similar resolution in 
the magnetic field calculations. For modelling satellite observations, this 
provides a distinct advantage over previous systems which had resolution of 2°-3° 
In latitude and 15° In longitude. 

At present, because of numerical quadrature restrictions, the grid 
cell system does not extend beyond 85°N (centered dipole). An additional system 
for modelling the perturbation magnetic field due to currents above 85°N is 
under development. This system replaces ’’cells'' by transpolar current "strips". 
In a manner similar to the cell system, aB due to current in each strip will be 
computed and stored (kernels). Superposition of a number of these strips, at 
different orientations, and weighted appropriately, will generate resultant polar 
cap current systems of sufficient complexity to satisfy most modelling needs 
(Figure 3). 


As in the case of the cell system, all quadrature associated with 
the solution of the Bio-Savart law is carried out only once for a basic strip 
system. All other systems can be generated from this set of kerr«ls. 

3.3 Data Analysis 

a) A study of the orientation of the polar cap residual vectors in the H’D’ 
plane has b^^en undertaken. Polar transits satisfying high IMF stability 
criteria, and having "well-ordered" vectors have been examined. To date, 
"well-ordered" remains subjective, and In^lles perturbation vectors which 
are reasonably uniform in magnitude and direction, directed sunward in 
the northern hemisphere and anti sunward In the southern hemisphere. Only 
data from transpolar passes whose polar cap vectors extend at least 5° on 
either side of the noon-midnight meridian, and with a total extent of at 
least 15° poleward of the auroral oval are considered. Figure 4 shows 
examples of acceptable transits for both hemispheres. Most of these 
passes occur at times when the Z-con^orient of the IMF Is negative, and 
their frequency Increases with Kp. 

If the average orientation of polar cap vectors in this data subset is 
defined as the average angle of those vectors falling In the central 60 
seconds of data (I.e., ±30 s from the noon-nrtdnlght meridian), then a 
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correlation coefficient of 'v4).7 U obtained betifeen this average vector 
direction and the Y-component of the INF. This Is true In both hemi- 
spheres . 

A similar attempt to examine the average vector magnitude has been some- 
what Inconclusive* Indicating* to date* only that this quantity Is larger 
In the sunlit (southern) polar cap than In the dark polar cap. 

b) A statistical analysis of the fluctuation parameter (FP) (a measure of 
tte high frequency content of the magnetic (^rturbatlon vectors; see 
teport #1) has been performed using biro months of Nagsat data (November 
1979, January 1980). Data have been sorted Into bins 1^ of latitude 
by 1 hour of magr^tlc local time* and examined for the presence of 
significant FP as a function of Kp (3 bins: 0 to 1+; 2- to 3+, >3+) 
and IMF By (By < 0; By >0). The selection criteria were such that 26 
days are Included In the lew/ Kn case* and only 7 In the high Kp case. 

Figure 5 shows the frequency of occurrence plots for thg Kp analysis. 

Now that Nagsat data up to May 1980 have been distributed, this analysis 
will be performed on tSe expanded data set. In the meantime, no defini- 
tive Interpretation cf the results has been made. 

3.4 Cooperative Activities 

a) Discussions with other Nagsat Investigators brought to light a problem 
with the method of computing Nagsat magnetic field residuals with accept- 
able accuracy. At NRCC (as at a nunter of institutions), the model 
field was Initially confuted at one-minute intervals along the orbit, 
and then Interpolated from these values to 1 -second values* using a 
cubic spline interpolation scheme. Because of rapid changes in the field 
in the neighbourhood of the poles, this method introduces oscillations 
there up to 40 nT in the e and 41 con^onents. Use of a 5-point Lagrange 
interpolation code leads to even higher amplitude oscillations. However^ 
altering the scheme so that the model field Is calculated every 15 s, 
and then interpolated to 1 s values, removes the oscillations. Figure 6 
Is a comparison of the two methods of Interpolation. 

b) A cooperative study with Rothwell et al . (University of Southampton) has 
been undertaken. Auroral video data from Longyearbyen, Spitzbergen for 
January 15, 1980 at 0715 UT indicates an arc near noon MLT which Is 
coincident with a Nagsat orbit. The analysis is incon^lete at this time. 

(See Daytime auroral motions imaged with an all -sky low light level TV 

camera In Spitzbergen, P.A. Rothwell, B. Lanchester, R.W. Thomas, lAGA* 1 

Edinburgh, 1981, paper 2C.08). 

i 

3.5 Data Processing Status 

We have at present processed 38 days of Nagsat data from the interval 
launch through January 1980. We are ready to process the remainder of the Nagsat 
data. 

1 
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4 . Significant Results 

The unprecedented precision of the vector data set fr«n Magsat has 
made It necessary to undertake an extended exploratory phase for data analysis 
procedures, modelling techniques and phenomenology. Although a number of signi- 
ficant preliminary results have been obtained, as outlined above, the studies 
are still In progress and have not yet been submitted for publication. 

5. Publications 

a) Comparison of quiet day B-fleld variations In the summer and winter auroral 
zones . 

M.O. Wilson, T.J. Hughes, J.R. Burrows 
Spring AGU 1981, paper (3^45. 

b) The rc e of crustal magnetic anomalies In modelling Ionospheric current 
systems. 

T.J. Hughes, D.D. Wallis 
Spring AGU 1981, paper GP44. 

Papers are in preparation for the Magsat special issue of Geophysical 
Research Letters. 

6 . Problems 


We reiterate our concern expressed in Report #1 regarding the high 
frequency noise that has been introduced into CHRONFIN data by the fine attitude 
calculation. We had hoped to study currents at the maximum spatial resolution 
permitted by the Magsat data rate. Introduction of this ''solution noise" has 
reduced this resolution considerably. 

In an attempt to understand the source of residual errors in the fine 
attitude solutions and perhaps develop a 'work-around* solution to them, we have 
made a limV.ed study of the attitude flags and corresponding attitude solutions. 
There appears, on occasion, to be an unexpected Inconsistency in them. An 
example is described in Appendix A. This situation has inhibited our attempts 
to develop an algorithm for correcting motion model segments of the fine atti- 
tude solution. The study is temporarily suspended. 

7. Recommendations 

No further recommendations at this time. 

8. Conclusions 


Although there have been no definitive results to date, preliminary 
investigations are encouraging. Analysis should proceed reasonably quickly now 
that data from most of Magsat' s lifetime has been distributed. 
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Figure Captions 

Fig. 1. An example of the NRCC standard UT plot of magnetic perturbations In 
the H'D'U* coordinate system. 

Fig. 2. An exan^le of the NRCC standard polar plot of magnetic perturbation 
vectors . 

Fig. 3. Schematic of the strip current system for polar cap currents. In prac- 
tice* up to 45 strips and 45 rotations may be used. 

Fig. 4. Examples of “well-ordered" polar cap perturbation vectors. 

Fig. 5. Polar plot of frequency of occurrence of significant fluctuation para- 
meter (FP > 8 nT) for combined north and south hemispheres, for November 
1979 and January 1980. Contours are drawn at 20%, 40%, 60% and 80%, 
with the 60% contour heavier. No coverage occurs where the plot Is 
blank . 

a) 0 < Kn i 1+ 

b) 1+ < Kp < 3+ 

c) Kp > 3?. 

Fig. 6(a). Magnetic perturbations relative to the M6ST 06/80 model field calcu- 
lated every minute along the orbit and interpolated to values 1 second 
apart. Note the H' and O' conponents between minutes 46 and 48. 

(b). The perturbations relative to the model field conputed every 15 seconds 
and then Interpolated to 1 s values. The character of all 3 compo- 
nents has been changed significantly from that In Figure 6(a). 

Fig. 7. Plot of magnetic perturbations for day 316 1979 near 0440 UT. See 
Appendix A for discussion. 


Appendix A 

There appears* on occasion, to be a problem with the attitude solutions 
and their flags. For example, day 316 (1979), near 04:40 UT (Figure 7): there Is 
a distinct neget1ve«go1ng step of 95 nT In the X-component, and 11 s later, a 
return (positive- going) step of 149 nT. This latter step Is flagged as a shift 
out of the motion model (flags 7608 to 1022), whereas the first shift occurs 
In the middle of a period of motion model. This Is not expected from our under- 
standing of the attitude solutions and the flags. That Is, we do not expect 
to see such shifts In the data except when there has been a change In type of 
attitude sol'.tlon. It Is not known at this time If other such unflagged shifts 
aie in the data. 
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